Impurity effects on s+g-wave superconductivity in borocarbides Y(Lu)Ni 2 B 2 C 
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Recently a hybrid s+g-wave pairing is proposed to describe the experimental observation for a nodal 
structure of the superconducting gap in borocarbide YM2B2C and possibly LUN12B2C. In this paper 
the impurity effects on the s+g-wave superconductivity are studied in both Born and unitarity limit. 
The quasiparticle density of states and thermodynamics are calculated. It is found that the nodal 
excitations in the clean system are immediately prohibited by impurity scattering and a finite energy 
gap increases quickly with the impurity scattering rate. This leads to an activated behavior in the 
temperature dependence of the specific heat. Qualitative agreement with the experimental results 
is shown. Comparison with d-wave and some anisotropic s-wave studied previously is also made. 



PACS numbers: 74.20.Rp, 74.62. Dh, 74.70.Dd 
I. INTRODUCTION 

The superconductivity in rare earth nickel borocar- 
bides RNi 2 B 2 C (R=Y, Lu, Tm, Er, Ho, and Dy) is of 
great interest in recent years. 1 ' 2 Among them, the two 
nonmagnetic borocarbides YM2B2C and LUM2B2C have 
relatively high superconducting (SC) transition temper- 
atures, which are 15. 4K and 16. 5K, respectively. Al- 
though they were initially understood by an isotropic 
s-wave pairing, 3 recent various experimental results, in- 
cluding specific heat, 4-7 thermal conductivity, 8,9 Ra- 
man scattering, 10,11 NMR relaxation rate, 12 photoemis- 
sion spectroscopy, 13 scanning tunneling microscopy and 
spectroscopy, 14 have shown that they may be another 
class of nodal superconductors. For example, both the 
VW dependence of the specific heat and H linear de- 
pendence of the thermal conductivity in the vortex state 
indicate the existence of nodal excitations. 4-6 In partic- 
ular, compelling evidences are presented by Izawa et al. s 
based on the angular dependent thermal conductivity 
measurements in a magnetic field that the gap function 
of YN12B2C is highly anisotropic with anisotropy ratio 
less than 10 -2 , i.e., it has essentially nodes and more- 
over, the nodes are point-like which are located along 
[1,0,0] and [0,1,0] directions. Due to Doppler shift in the 
presence of supercurrent flow, 15 the zero energy quasi- 
particle density of states (DOS) becomes finite and is 
sensitive to the relative orientation between the magnetic 
field and nodal direction. Thus the nodes can be unam- 
biguously extracted from the angular modulation of the 
thermal conductivity. Based on the same physical mech- 
anism the angular dependent specific heat in YNi2B2C 
was measured by Park et al. 7 and the conclusion for the 
existence of the point nodes can be also drawn. In ad- 
dition, the thermal conductivity in LUM2B2C was mea- 
sured by Boaknin et al. as a function of temperature 
and field strength. 9 They claimed that the gap minimum 
is at least 10 times smaller than the gap maximum and 
possibly goes to zero at nodes. Thus YM2B2C (and pre- 



sumably LuNi2B2C) becomes the first superconductor in 
which the gap function with point nodes is clearly iden- 
tified. 

Theoretically, Maki et al. have proposed that the so 
called s+g-wave spin singlet gap function, i.e., 16 ' 17,8 



A k = A(l - sin 4 6>cos40)/2 = A/(k) 



(1) 



can describe the experimental ob- 
servation in Y(Lu)Ni2B2C. Here 9, <fi are the polar and 
azimuthal angles of k, respectively. The second 'g-wave' 
component is given by a fourth degree Ai s basis function: 
fc 4 + ky — 6fc 2 fc 2 in tetragonal crystal symmetry. In the 
gap function the amplitudes of s- and g-components are 
assumed to be equal. Thus 4 (and only 4) point nodes 
at = ir/2 and <j> — 0, tt/2, it, 3ir/2 are realized, which is 
exactly what has been observed experimentally. 8 

So far, a microscopic theory for the s+g-wave pairing 
is not available which might be complicated due to the 
complex Fermi surface of borocarbides 18 and the possi- 
bility of strongly anisotropic Coulomb interactions. A 
phenomenological theory has been constructed by two of 
us, 19 which shows that a stable coexistence of s- and g- 
wave with the fine tuning, i.e., the equal amplitudes of 
them are possible. On the other hand, the correctness 
of the s+g model could be checked by comparing the 
theoretical predictions with the available experimental 
data. Up to now, quite a few physical quantities based 
on the s+g model such as thermal conductivity, 8 Raman 
spectra, 20 sound attenuation 21 etc. have been calculated 
and good agreement with the experimental results is ob- 
tained. This in turn stimulates us to continue studying 
more properties. 

In this work the impurity effects on the s+g 
model will be investigated. Experimentally there have 
been extensive studies on Pt-substituted borocarbide 
Y(Nii_ a; Pt :c )2B2C. The ultrahigh-resolution photoemis- 
sion spectroscopy was performed by Yokoya et al. 13 They 
found a small but significant difference for the spectra in 
the SC states between x = and x = 0.2 compound, 
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which is contributed to the change in anisotropy of the 
SC gap: a strong anisotropic gap in x — and an almost 
isotropic one in x = 0.2. That is, the gap anisotropy 
is wiped out by introducing impurities. Consistent re- 
sult can be obtained from the electronic specific heat 
G measurements by Nohara et al. 5 The temperature- 
dependence of G (under zero field) shows a power law 
behavior for x — but a thermally activated behavior 
for x = 0.2. Very recently the angular dependent ther- 
mal conductivity for Y^ii-^Pt^I^C with x — 0.05 in 
the SC state was measured. 22 There is no angular modu- 
lation left. This shows that the gap anisotropy is nearly 
completely suppressed by only 5% Pt substitution. It is 
the purpose of this work to see how soon a full energy 
gap will be opened by the impurities for the s+g-model 
through the calculation of the DOS. Furthermore, the 
thermodynamics like specific heat will be calculated. The 
crossover from a power law behavior in a pure system to 
an activated behavior in the presence of impurities will 
be exhibited. 



II. FORMALISM 



Above & n = uj n + zSo, Ak = Ak -I- Si, and (• • •) means 
angular average J dO/(47r). It deserves to point out that 
the off-diagonal self-energy Si is zero for d-wave but fi- 
nite for the current s+g-wave. This difference will lead 
to different behavior of their respective DOS functions. 

The gap function is connected to the matrix clement 
gi2 by A k = TJ2 iuJn J2w Mck'#i2(k', iu> n ) with Vkk': the 
pair potential. It is often assumed Vkk' = — V/(k)/(k'), 
which is valid for \u) n \ < ujd (a cut-off). Then the follow- 
ing equation for the gap amplitude A is obtained: 



A = 2ttTV 




/(k)(A k + £i) 



v /K + ^ ) 2 + (A k + S 1 )2 



(5) 



where V — VNq is a dimcnsionless interaction constant. 

By solving Eqs. (3)-(5) self-consistently, we can obtain 
the renormalized A with impurity concentrations at all 
temperatures and further calculate the observable quan- 
tities like specific heat. In the following we adopt the 
impurity scattering rate parameter in the normal state 



r = r„/(i 



which reduces to T u /c in Born limit 



(c 3> 1) and r„ in unitarity limit (c = 0). 



The effects of nonmagnetic impurities are treated by 
the self-consistent T- matrix approximation. 23 For benefit 
of readers, we briefly summarize the formulation in the 
following. To begin with, the Matsubara Green's func- 
tion in Nambu representation is written as (quantities 
with a hat represent matrices) 



<?(k, iuj n ) = (iu n a - £k03 - A k <ri - E)~ 



(2) 



where ui n = (2n+ 1)ttT (T: temperature) is the Matsub- 
ara frequency, £k the electron energy measured from the 
chemical potential, Ak the gap function which has been 
assumed to be real, ao and 0-1,3 are unit and Pauli ma- 
trices, respectively. The self-energy E(iw„) = T u T(iu n ) 
with r„ = Ni/(irNo). Here Ni is the total number of 
impurities and N is the zero energy DOS per spin for 
the spectrum £k- The T-matrix T satisfies T(iu> n ) = K + 
KG(iu n )T(iuj n ), where G{iu n ) = (1/tvN ) J2k 5( k , iw n ). 
K is the T-matrix in the normal state: K = —a^/c with 
c = cot#o- And #o is the s-wave phase shift controlled by 
the strength of a single impurity potential. 24 

The matrix A (A = E, T, G) is decomposed into A = 
Aq<7q + A\<j\ + Then one obtains G3 = if as- 

suming particle-hole symmetry. 23 Further T3 and corre- 
spondingly S 3 vanish if we only concern Born (c » 1) 
and unitarity (c = 0) limit. The left So,i are obtained as 
follows 



Go 



c 2 -G 2 +G 2 



2 ' 



where 



Go 



-lUJr, 



4^ 



Si — r„ 



Gi 



-Gi 



c 2 - Gq + G\ 



-A, 



2 ' 



+ A 2 



4^ 



+ A 2 



(3) 



(4) 



III. DENSITY OF STATES 

We first consider the DOS, which is given by (per spin) 

N(E)/N = -TmG (iu n )\ iUn ^ E+m . (6) 

It is noticed that if one scales the input parameters F and 
E by A, the DOS can be obtained by only solving Eqs. (3) 
and (4) without explicit calculation of A for the moment. 
The results for a few ratios T/A are shown by Fig. 1 
in both Born and unitarity limit. (Note that the scaled 
r/A is still a monotonous function of T.) From the figure 
a few features are clearly seen: i. There is no essential 
difference between the results in both limits, as seen by an 
example in Fig. Ib. Numerically we have — Gq + G\ ~ I 
in unitarity limit. Thus Eq. (3) can be rewritten into 
a unified form: E = TGo, Ei = — YG\ in both limits, 
leading to nearly the same DOS. ii. An energy gap opens 
up immediately when r ^ and grows quickly with T. 
As seen from the inset of Fig. la, a relatively big energy 
gap uj g <~ T is already present for small r/A = 0.02. 
Actually uj g is well approximated by u) g = T/(l + 2r/A). 
For intermediate r/A = I, the gap anisotropy is nearly 
suppressed and the DOS approaches the familiar shape 
of an isotropic s-wave gap. All the above results are quite 
different from those for d-wave. 25-31 For example, a finite 
DOS is induced by impurities around zero energy for d- 
wave. Also, a resonance behavior at low energy exists in 
unitarity limit for d-wave, which does not appear in the 
s+g-state. We point out that the main result seen from 
the DOS based on the s+g model, i.e., the gap anisotropy 
is suppressed by impurity scattering, is consistent with 
the conclusion drawn from the photoemission data. 13 
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FIG. 1. The quasi-particle DOS for s+g-wave. (a) The 
results for various T in Born limit. The low energy region 
is enlarged in the inset, (b) Comparison between unitarity 
and Born limit for a fixed T/A = 0.5. Only subtle differences 
exist in both limits. 

At this stage, we would supplement a remark. Previ- 
ously, in order to clarify the role of sign change in a gap 
function, some hypothetical anisotropic s-wave gap func- 
tions e.g. | cos 2(f>\ , 26 which vanishes at the same nodes as 
d-wave but never change sign, have been studied in the 
presence of impurities. The obtained DOS for r ^ has 
qualitatively similar properties to that shown here. This 
is understandable in the sense that the current s+g-wave 
(without sign change) is essentially an anisotropic s-wave. 
On the other hand, however, we would point out an im- 
portant difference between the current s+g-wave gap and 
the one ~ |cos2(/>|: The former contributes only a few 
point nodes, while the latter gives line nodes, on a three- 
dimensional Fermi surface. This will lead to notable dif- 
ferences on some observable quantities, e.g., the angular 
dependent thermal conductivity, which has actually been 
used to identify whether the SC gap in YN12B2C has re- 
ally point or line nodes. 8 As far as the DOS is concerned, 
a quantitative difference is that the opened gap ui g in- 
creases with T for s+g-wave much faster than that for 



| cos 2<p\ since point nodes are easily removed by impu- 
rity scattering. This is consistent with the experimental 
observation that YNi 2 B 2 C in the SC state is very sensi- 
tive to impurities as mentioned above. 



IV. TRANSITION TEMPERATURE, SC ORDER 
PARAMETER AND THERMODYNAMICS 

In this section we focus on the thorough solution of 
Eqs. (3)-(5). Considering the limit A — > 0, we first 
extract the transition temperature T c analytically. In 
both Born and unitarity limit, we obtain the following 
equation which is different from the standard Abrikosov- 
Gor'kov (AG) one 32 



-ln^ = 



(J 2 (k)) - (J(k)) 2 

</ 2 (k)) 



, (7) 



where T c q is the transition temperature at T = 0, and 
ip(z) is the digamma function. The prefactor p — 

^ ^ k (f2(k)) k ^ — 0-169 for the current s+g-wave. Note 
that for the usual nodal superconductors with only a sin- 
gle angular momentum pairing like d-wave or pure g- 
wave, the prefactor is unity, i.e., Eq. (7) becomes the 
standard AG equation. 27 ' 28 Therefore, for d-wave there 
exists a critical T: T c = nT c0 /{2e^) ~ 0.882T c0 at which 
T c vanishes, whereas for s+g-wave T c = oo. In the limit 
T -> 0, Eq. (7) gives 



Tc = T c0 - -pT . 



(8) 




r/A, 



00 



FIG. 2. T c /TrM (solid lines) and A(F, T = 0)/A 00 (dotted 
lines) vs. r/A o for s+g- and d-wave. Note that for each 
quantity the results in two limits are the same for s+g-wave. 
For d-wave the quantity A(r,T = 0)/Aoo is shown in unitar- 
ity limit. 28 
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The dependence of T c on T is shown in Fig. 2, where 
r is scaled by A 00 = A(T = 0,T = 0). The result for 
d-wave is also shown for comparision. From Eq. (8) we 
have the linear relation at T — > 0: T c /T c0 = 1 — qT / A o 
with q — 0.366 for s+g-wave and 1.68 for d-wave. Here 
we have used the BCS ratios Aoo/T c o ~ 2.76 and 2.14 
for s+g- and d-wave, respectively. 

We try to evaluate T from the experimental data on 
T c . For Y(Ni 1 _ x Pt x ) 2 B 2 C, T c = 12. IK at x = 0.2, 13 
i.e., T c /T c0 ~ 0.786. Then we have r/A 00 ~ 0.94. For 
x = 0.05, T c = 13. IK, 22 corresponding to r/A 00 ~ 0.55. 

The SC order parameter A(r, T) is solved numerically. 
In what follows we give the results only in Born limit. 
Those in unitarity limit are practically the same. In 
Fig. 3, A(r,T) is given as a function of T for a few T. 
In addition, A(r,T = 0)/A o is shown in Fig. 2, which 
decreases slowly (slower than T c /T c0 ) with T. 




FIG. 3. A(r,T)/A o as a function of T/T cQ at different 
F/Aoo for s+g-wave. 
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FIG. 4. The entropy S/S„ as a function of T/T c o for a few 
r/Aoo- The inset shows S/(jT c0 ) vs. T/T c0 . 



We further calculate the thermodynamics, 
tropy S is given by 



The en- 



/OO 
dE N(E) {f(E) Hf(E)} 
-00 

+ [1 -/(£)] ln[l -/(£)]} (9) 



with Fermi function ,f(E) = 1/(1 
have shown the entropy S/S n vs. 



- e^ E ). In Fig. 4 we 
T/T c q for a few T, 
where S n = 7T with 7 = ^tt 2 N is the entropy in the 
normal state. In the inset the same figure is reploted 
but S is scaled by a constant 7X^0. This is helpful to 
understand the specific heat shown below. 

The specific heat is given by C = Tff , which can be 
explicitly expressed as 



C_ 



2tt 2 T 3 



+ 



_6_ r 



dE 



N 
N(E) 



dT 



N 



ln(l + el m ) - 



f3E 



(10) 



where C„ = 7T is the specific heat in the normal state. 
In Eq. (10) the second term is contributed by the T- 
dependence of the DOS. It is important only at high 
temperatures because at low T, about T < 0.2T c o (see 
Fig. 3), A(r,T) is nearly T-independent and also the 
DOS. For an instructive look at how the DOS changes 
with T in the whole range, we have shown in Fig. 5 an 
example for r/A o = 0.2. 




E/A QQ 



FIG. 5. The DOS with change of T for r/A 00 = 0.2 (with 
corresponding T c /T c o = 0.934). It is seen that the DOS at 
T/T c0 = 0.3 is only slightly changed from that at T = 0. With 
increasing T, the bump at relatively high energy dissolves 
gradually, and the peak becomes narrow and simultaneously 
shifts towards zero. When T — » T c , the DOS is close to an 
inverse square root function valid for an isotropic s-wave gap. 
Finally it will become the straight line "1" at T — T c . 
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FIG. 6. The specific heat C/C„ vs. T/T c0 for 
T/Aoo = 0, 0.2 and 1.0. The low temperature region is en- 
larged in the inset. 

The specific heat C/C n vs. T/T c0 is shown in Fig. 6. 
We mainly concern the low temperature behaviors, which 
are enlarged in the inset. The notable observation is that 
the T-dependence changes from a power law behavior in 
the pure case into an activated behavior in the presence 
of impurities. This is due to the opening of a full gap by 
impurities as seen from the DOS. Thus the low energy 
quasi-particle excitations are exponentially small. 

The above point becomes more transparent if we 
present the analytical results for the specific heat at low 
T where only the first integration in Eq. (10) contributes. 
In the pure case, the DOS is linear at low energy E <C A: 

Tvp = f f ■ For about °- 1A / T > 5 ' Le -> T < 0-05T c0 , 
this linear formula can be used throughout the whole 
range in the integration. Then we have 



an 



i.e., a linear T-dependence (or T 2 dependence for C). 
When r > 0, the lower limit of the integration becomes 
u g . In view of sech 2 (£/2T) ~ 4exp(-£y"T) at low T 
(about ujg/T > 5), an exponential factor ~ cxp(— u> g /T) 
is naturally expected after the integration irrespective of 
the concrete form of N(E). Explicitly, for T < T < A 00 
we have obtained (u> g ~ T) 



C_ 

C„, 



3T 5/2 



^3/2 T 3/2 A 

277 \ 
X(1 + T )e 



ln(l 



A, 



r/T 



(12) 



Experimentally the specific heat was measured by No- 
hara et al. 5 for Y(Nii_ x Pt x ) 2 B 2 C with x = and 0.2. 
By fitting the experimental data they have exhibited a 



power law relation C cx T n with n ~ 3 for x — com- 
pound and C cx exp(— const. /T) for x = 0.2. The essen- 
tial features of the experimental data in both pure and 
impure systems have been reproduced by our theoretical 
results, although some differences exist for the accurate 
T-dcpcndcnces between theory and experiment. 

In addition, the specific heat jump at the transition 
temperature can be analytically derived. For a pure sys- 
tem, we have 



C — Cn 



1.426 



T=T C( 



(/ 2 (k)) 2 
</ 4 (k)) 



0.885 . 



(13) 



For T > 0, the formula is algebraically more complicated: 



C — C n 



C n 

7C(3) 



= 12(1 + a) 2 



T=T, 
b 



1 - 



— X^(- + X) 



TV 

T 



^ 1 \-+x) 



ax 



^+^ 1 H\+x)-2x-^( 1 -+x)-^{ 1 -)^ 



(14) 



where x = T/(2irT c ), a = (w 2 ), b = (w 4 ) with w = 
sin 4 6 cos 40, and tp^(z) (i = 1, 2, 3) are the poly-gamma 

C-C n 



functions. From Eq. (14), we obtain — ^j- 

0.901(1.016) for T/Aoo = 0.2(1), which are consistent 
with the numerical values in Fig. 6. 



V. CONCLUSION 

In the present paper we investigate the impurity effects 
on s+g-wave superconductivity in YNi 2 B 2 C and possibly 
LuNi2B2C. Unlike in the usual nodal superconductors 
we have discovered that the results in Born and unitar- 
ity limit are practically the same. There is no resonant 
scattering associated with nonmagnetic impurities. Also 
unlike in the usual nodal superconductors, the quasi- 
particle energy gap opens up immediately by impurity 
scattering and grows quickly with the scattering rate. In 
other words the nodal excitations are prohibited imme- 
diately when impurities are introduced, which is clearly 
seen from the quasi-particle density of states. Indeed we 
could evaluate that the energy gap opens up even for 
Y(Nii_ ;E Pt a; )2B2C with x = 0.01. From the value for T 
at x = 0.05, we estimate T/Aoo ~ 0.1 for the crystal with 
x = 0.01, which is substantial for a gap opening. There- 
fore the s+g-wave superconductivity is very sensitive to 
impurity scattering. 

Also we have calculated the specific heat, which shows 
a change from a power law dependence of T in a pure sys- 
tem into an activated behavior in the presence of impu- 
rities. Compared with the experimental data by Nohara 
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et al., 5 the theoretical results have reproduced the main 
features. The accurate comparison needs improvements 
in both theoretical and experimental side. 

Finally we point out that the sensitivity to impurities 
predicted by s+g-model is confirmed by the angular de- 
pendent thermal conductivity measurements on the sin- 
gle crystal Y(Nii_ x Pt x ) 2 B 2 C with x = 0.05, 22 where the 
gap anisotropy has been cleared out. In a forthcoming 
paper we shall discuss other properties such as supcrfluid 
density etc. within the same formal treatment. 
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